1 



SPECTRAL ANALYSIS IN THE VIRTUAL OBSERVATORY 



Thomas Rauch 1 and Iliya Nickelt 2 

'institute for Astronomy and Astrophysics, Kepler Center for Astro and Particle Physics, Eberhard Karls University, 

Sand 1, 72076 Tubingen, Germany 
2 Astrophysical Institute Potsdam, Department of E-Science, An der Sternwarte 16, 14482 Potsdam, Germany 



ABSTRACT 



In a collaboration of the German Astrophysical Virtual 
Observatory (GAVO) and AstroGrid-D, the German As- 
tronomy Community Grid (GACG), we provide a VO 
service for the access and the calculation of stellar syn- 
thetic energy distributions (SEDs) based on static as well 
as expanding non-LTE model atmospheres. 

At three levels, a VO user may directly compare ob- 
served and theoretical SEDs: The easiest and fastest way 
is to use pre-calculated SEDs from the GAVO database. 
For individual objects, grids of model atmospheres and 
SEDs can be calculated on the compute resources of 
AstroGrid-D within reasonable wallclock time. Experi- 
enced VO users may even create own atomic-data files 
for a more detailed analyses. 
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1. INTRODUCTION 



Spectral analysis of hot, compact stars by means of 
Non-LTE (local thermodynamic equilibrium) model- 
atmosphere techniques has for a long time been regarded 
as a domain of specialists. In contrast to the assump- 
tion of LTE, where occupation numbers of atomic levels 
are determined by Saha equation and Boltzmann statis- 
tics, radiative and collision al transitions h ave to be con- 
sidered in detail (see, e.g., IWernerl Il986l) . This makes 
the calculation of an elaborated model and its spectral 
energy distribution (SED) time-consuming, although it 
can be performed on presently available PCs. Faster and 
easier is the use of blackbody SEDs, tempting, e.g., old- 
fashioned users of photoionization codes to use such for 
the (bad) representation of the exciting source of an ion- 
ized gaseous nebula - spoiling the advantages of the best 
of these codes. 

Within the last four decades, Non-LTE model- 
atmospheres arrived at a high level of sophistication: 



fully line-blanketed models consider opacities of al l 
elements from H to Ni (lRauchll2003llRauch et al.ll2007l) . 
Within the GAVCQ, we provide access to SEDs of such 
models (Sect. [2j. 

In general, spectral analysis requires grids of SEDs over 
a wide parameter range (effective temperature T e g, sur- 
face gravity g, element abundances). Since their calcula- 
tion may exceed the available capabilities, model atmo- 
sphere calculations are performed on compute resources 
of AstroGrid-D. 



2. THEOSSA - MODEL SED ON DEMAND 



We have created the VO service TheoSSA (Theoretical 
Simple Spectra Accesfl At three levels it provides 
access to SEDs that are in a pilot phase calculated 
by the Tubin gen Mo del Atmosphere Packag e TMA/Q 
(IWerner et all 120031: iRauch & Deetjenl l2003h for hot, 
compact stars only: 



Fast and Easy 

An unexperienced VO user (no detailed knowledge about 
the model -atmosphere code necessary) may use pre- 
calculated SEDs, interpolate in between. 

Individual 

An interested VO user can calculate own SEDs for the 
analysis of special objects. 

Experienced 

The VO user can define own model atoms and upload 
them for the comparison with other codes, etc. 

Figure Q] shows the scheme of the dataflow within 
TheoSSA. If individual parameters are requested by a 
VO user, the GAVO database is checked for suitable 
SEDs (parameters match within error limits). If there is 

1 http://www.g-vo.org 

2 http://vo. ari.uni-heidelberg.de/ssatr-0.01/TrSpectra.jsp 
3 http://as tro.uni-tuebingen.de/~rauch/TMAP/TMAP.html 
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TheoSSA request: T e s, log g, X[ 



check GAVO database: 
requested parameters matched 
within tolerance limits? 



yes: offer existing SED 
I 



3_ 



no: calculate new model and SEDs 
and ingest into GAVO database 



3L 



request exact parameters 




accept, retrieve SED 



GAVO database @ AR1 

Teff, log<7,*i 
metadata 
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T cff , l0g<?,Xi 

models 
atomic data 
frequency grids 
^ SEDs j 



Figure 1. Data flow within TheoSSA. ARI: Astronomisches Recheninstitut (Heidelberg, Germany), IAAT: Institut fur 
Astronomie und Astrophysik Tubingen (Germany), X\: mass fraction of species i E [H, He, C, N, O]. 



no match, the VO user is guided to With this 

WWW interface, the VO user may calculate an individ- 
ual model atmosphere, requesting effective temperature, 
surface gravity, and mass fractions of H, He, C, N, and 
O (more species will be included in the future). For 
this calculation, standard model atoms are used which 
are provided within the Tubingen Model-Atom Database 
TMAD^. Since the VO user can do this without detailed 
knowledge of the programme code working in the back- 
ground, the access to individually calculated SEDs is as 
simple as the use of pre-calculated SEDs - however, the 
calculation needs some time (depending on the number of 
species considered, the wall-clock time is ranging from 
hours to a few days). Standard SEDs of all calculated 
model atmospheres are automatically ingested into the 
GAVO data base and, thus, it is growing in time. 

If a detailed spectral analysis is performed, an experi- 
enced VO user may create an own atomic data file tai- 
lored for a specific purpose considering all necessary 
species and calculate own model atmospheres and SEDs. 



3. THEOSSA @ GRID 

AstroGrid-D^ is a German research and development 
project to build an astronomical research infrastructure. 

4 http://as tro.uni-tuebingen.de/~TMAW/TMAW.shtml 

5 http://astro.uni-tuebingen.de/~rauch/TMAD/TMAD.html 

6 http://www.gac-grid.de/ 



Using the Globus Toolkit (GT4) middleware, AstroGrid- 
D embeds existing computational facilities, dedicated re- 
sources and specialized hardware, such as robotic tele- 
scopes. Almost a hundred users from the German astro- 
nomical community use its compute and storage facili- 
ties. The project also resulted in several scientific appli- 
cations and novel Grid services, such as the information 
service Stellari^\ 

The TMAP model calculations are carried out on the re- 
sources of AstroGrid-D to ensure a reliable service, even 
for a high number of requests. For job submission, a 
script package is used that allows to submit jobs to com- 
putational resources automatically. This script package 
was originally developed for a different AstroGrid-D use 
casfl but can be easily adopted to other implementations 
of so-called atomic grid jobs. The term refers to a spe- 
cific type of grid application, where all input data is trans- 
ferred together with the software without need for inter- 
process communication. The grid is thus used as a task 
farming mechanism. But even if this is only a compara- 
tively simple application for the use of the flexible GT4 
middleware, we found task farming to be an important 
requirement of many use cases, such as TMAP. 

The transfer of the TMAP model software and the input 
data, as well as the transmission of the results is part of 
the standard job submission process in Globus Toolkit. 
We make use of the GT4 web services and control the 

7 http://stellaris. zib.de/ 

8 http://www.gac-grid.org/project-products/Applications/ 
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process by a template written in a JSDL (Job Submission 
Description Language^. Whenever a user requests new 
model data to be calculated, the template is adopted to 
the specific case and a target machine is selected from a 
given list. Upon arrival, the results are passed on to the 
TheoSSA service to notify the user. 



4. CONCLUSIONS AND FUTURE PLANS 



TheoSSA is already fully functional and provides easy 
access to SEDs of hot stars - use them! E.g. the use 
of blackbody SEDs to represent stars in photoionization 
models belongs to the last millennium and is not adequate 
anymore. 

The database of complete model atoms for other species 
(Ne, Na, Mg, Si, ...) up to the iron group will be ex- 
tended. TMAW can already consider opacities of all ele- 
ments from hydrogen to germanium. 

The addition of other model codes for other types of 
stars will in future further improve the benefit for VO 
users. This is an invitation and challenge for other stellar- 
atmosphere working groups to contribute. 

Precise spectral analysis requires extended grids of 
elaborated model atmospheres. In th e framework of 
GRID computing jFosteii I2005L 120061) . the calculation 
of model-atmosphere grids and flux tables, e.g. via 
TheoSSA, is an excellent application to efficiently calcu- 
late synthetic spectra on reasonable time scales. 

One approach is already in preparation, using the newly 
developed HotBlast Non-LTE code for spherically ex- 
panding stellar atmospheres (Koesterke priv. comm.). 
HotBlast uses as an input the atmospheric structure of 
the static TMAP model atmospheres to simulate the at- 
mosphere below the wind region. 
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